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GRILLAGE BEAMS ON ELASTIC 
FOUNDATION! 


Bernard Drapkin# 


SYNOPSIS 


The geometric and elastic quantities influencing the maximum moment 
under a concentrated load for a beam on an elastic, homogeneous, isotropic 
halfspace, are represented by parametric curves. A numerical determina- 
tion of the spring constant, k, is effected for any particular case by equating 
the maximum moment determined from the curves to that of a beam resting 
on a spring foundation. The spring theory is subsequently applied using this 
readily available k; this would otherwise require experimental determination 
each time any physical quantity is varied, 

It is found that for the usual beams employed in practice there is a slight 
reduction in the maximum moment from that of the infinitely long beam. For 
foundations of high Young’s modulus it is quantitatively determined that the 
maximum moment is less than for foundations of low Young’s modulus. The 
effect of increasing the flexural rigidity of the beam is to increase the 
maximum moment. 

A simple formula is presented for steel beams on rock or average con- 
crete loaded by a single concentration at the center. If the ratio of Young’s 
moduli of beam and foundation is 15, the maximum 


Mom, of Inertia 
moment = 2/3 Load \3 173 ict = width’ The coefficient 2/3 is dimensionless. 


The maximum moment for elastic-rock foundations determined by the pro- 
cedure of this paper, is of the order of one-half that obtained when the usual 
practical design assumption of uniform pressure is employed. The maximum 
foundation pressure is of the order of twice the uniform foundation pressure. 
The assumption of uniform pressure therefore gives an understressed beam 
and an overstressed foundation under the concentration. 


INTRODUCTION 


This paper presents a solution of the grillage beam problem which is 
compatible with usual assumptions of the theory of elasticity and theoretical 
soil mechanics. 

The simplifying assumptions used in actual current and past practice to 
analyze and design these important foundation members are: 


a. Structural Engr., New York, N. Y. 

1, This paper is based on a dissertation submitted to the Polytechnic Institute 
of Brooklyn in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the field of applied mechanics. 
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a) The concentrated column load imposed on the beam is assumed to 
result in uniform pressure distribution on the supporting medium regardless 
of the relative stiffness of the foundation and beam, and independent of geo- 
metrical quantities such as beam length, beam width, foundation width and 
foundation depth; 

b) The selection of the beam is made by keeping the pressure below a cer- 
tain specified value and computing the maximum moment in the beam based 
on assumption (a) and the elementary beam theory. 


Most of the existing literature dealing with the subject of beams on elastic 
foundation assumes a known predetermined modulus of foundation k, the use 
of which simplifies the solution of the linear fourth order differential equation: 


4 
Ep I “a + ky = p(x).* The basic problem is the infinite beam on a spring 


foundation of modulus k, subjected to a single concentrated load. From this 
solution, by means of the superposition principle, solutions for various load- 
ing conditions have been determined. The direct solution of the finite beam 

on a spring foundation of modulus k is unwieldy due to the boundary conditions; 
the superposition procedure used extensively by Hetenyi2 is more elegant. 

A Fourier integral solution of the infinitely long beam carrying a concen- 
trated load and resting on an elastic foundation has been derived by Biot.3 
Biot obtained values for the foundation modulus, k, of the spring theory by 
equating the maximum moment in the infinite beam resting on the elastic foun- 
dation to the maximum moment in the same infinite beam on a spring foun- 
dation obtained from the elementary theory. Using this equivalent value of k, 
one can determine the moment, shear, deflection and pressure distribution 
along the beam using the elementary spring theory. 


Problem 


The basic structural element of a foundation grillage is a steel beam in 
full contact with the supporting medium and transmitting a heavy concen- 
trated load. Investigation of beams of finite length resting on an elastic foun- 
dation will be made to determine the maximum moment in the beam. From 
this result the equivalent spring constant k is determined. The effects of the 
length, width, flexural rigidity of the beam and the modulus of elasticity and 
Poisson’s ratio of the foundation on the equivalent spring constant for beams 
of any length, are presented in this paper by means of parametric curves. 
Comparison will be made with the results obtained by means of the usual de- 
sign assumption of uniform foundation pressure under grillage beams. 


Assumptions of the Theory of Elasticity used in Paper 


a) The supporting medium is an elastic foundation, homogeneous, isotropic, 
semi-infinite in depth and extent, without body forces, and characterized by 
the elastic constants, Young’s modulus E and Poisson’s ratio). 

b) The elementary beam theory is applicable. The beam has flexural 
rigidity Epl and negligible flexibility across its width of flange 2b. 


* See Notation 

2. “Beams on Elastic Foundation,” by M. Hetenyi, U. of Michigan Press, 
Ann Arbor, Michigan, 1946, p. 9. 

3. “Bending of an Infinite Beam on an Elastic Foundation,” by M. A. Biot, 
J. Appl. Mech., Vol. 4, 1937, pp. A1l-AT7. 
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c) The elastic foundation is always in contact with the beam and can 
support tension as well as compression. Young’s modulus is assumed to be 
the same constant in tension and compression. 

d) There are negligible shearing stresses at the contact surface. 


The first assumption in (c) is unnecessary when the contact pressure is 
sufficiently large to keep the beam in contact with the foundation. That is, 
the beam under its own weight induces a certain amount of compression over 
the entire contact surface; this initial compression may be everywhere 
greater than any induced negative pressure due to the concentrated load. 

The use of the elementary beam theory, within its limitations, is valid for 
grillage beams as for beams supported in other ways. The depth of the 
beams considered must be small compared with the length, since the elemen- 
tary theory assumes the shearing strains to have negligible effect on the axial 
strains. Straight prismatic beams are assumed to have proportions that pre- 
clude failure by twisting, lateral collapse and local wrinkling. The deflection, 
slope, moment and shear are proportional to the concentrated load since the 
basic differential equation is linear. The superposition principle and recipro- 
cal theorem are therefore applicable. 


Solution 


Biot* showed that the moment at any point in an infinitely long beam 
carrying a concentrated load P and resting on an elastic foundation is: 


M(x) = Pe (1)* 


The superposition principle and Maxwell’s reciprocal theorem are used to 
determine the free end effects for beams of finite length. Unknown moments 
Mo and unknown shears Q, are applied to those points on the infinite beam 
corresponding to the end points of the finite beam. The magnitudes and direc- 
tions of My and Q, are then determined so that the moments and shears in the 
infinite beam become zero at those points. The moment at any point in the 
finite beam is then expressed by adding the moment due to the concentrated 
load, the moment due to the moment Mp and the moment due to the shear Q,. 

Differentiating under the integral sign with respect to x in (1), the shear at 
any point in the infinite beam due to a concentrated load at the origin is 
obtained: 


4. “Bending of an Infinite Beam on an Elastic Foundation,” by M. A. Biot, 


J. Appl. Mech., Vol. 4, p. 1A, 1937. 
* See Notation 
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= ad (x 
V(x) Pc (x) (2) 


V dx V+dV 


Beam Element 


The beam sign convention used is: Shear taken positive when up on the 
left end of OB in Fig. 1. 


Fig. 2(a) is a limiting case of Fig. 2(b). 
M, = Pe = constant couple, as e~0, Pe-M,. 


Pe ~ Pof (x40) = Pe [F(x) - Fx + 


M(x) = wpe -Z(x +e) /, and in the limit as 
e 
=- 42) 
° 


M, at 0) 


Vix) = diz) 


(due 
M, at 0) 


M' and V' are the moment and shear at A due to P. 
M' = Pe 
We 


M" and V" are the moment and shear at A due to Mp1, Mog, Q,1» Qo2- 


(£10) +30) -% | + 


Now, M' + M"=0 and V' + V"'=0, since the moment and shear are 
zero at A. 


771-4 


Mag Mor nf 
AZLe*eZLE | | 
z 
t 
Fige 1 
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Fig. 2 
a) 
(5) 
(6) 
(7) 
(8) 
\ 


Pe J (f/2) + Qe - M/2=0 (9) 


Pe + Se 42) - %/2 + Moe - (10) 


Now, using J = x/c as dimensionless coordinate, 
Hence (9) and (10) velit 


PZ + #00) + - 1/2]=0 
c 


PF + - 1/2| + u, - =0 (10a) 
Cc 


This pair of simultaneous equations can be solved for Q, and M,: 


- 


Here the prime ' denotes differentiation with respect toy, 
i. e., 


An alternate form for My which is more suitable numerically when dealing 
with higher values of the argument {/2c is found by eliminating Q, between 


(9a) and (10a), giving: 


(13) 


The first term is the moment in an infinitely long beam and the remaining 
terms represent the effect of the shears and moments which have been de- 
termined by (11) and (12). The sum of these terms gives the moment at any 
station in the finite beam. It is necessary to evaluate the following types of 


integrals: 


Q=P. (2c) | Bc - 20)| g"(0) g"(Fe (12) 
beam 
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Integrals (14) and (15) are convergent whereas (16) is a divergent integral. 

From physical considerations @(7 ) represents the dimensionless quantity 
by which Pc is multiplied to obtain the moment at any point in an infinite 
beam due to a concentrated load at the origin. 

?' (J) represents the dimensionless quantity by which P is multiplied to 
obtain the shear at any point in an infinite beam due to a concentrated load, 
P, at the origin. @' ({) also represents the dimensionless quantity by which 
Mo, is multiplied to obtain the moment at any point in an infinite beam due to 
Mo at the origin. 

@'(y) represents the dimensionless quantity by which M,/c is multiplied 
to obtain the shear at any point in an infinite beam due to M, at the origin. 

In all numerical work the boundary conditions that the moment and shear 
are zero at the free ends of the beam, s= $/2c, are applied to give an addi- 
tional check on the work. 

The maximum moment, M, occurs under the load at J = 0 and is given by 
the expression: 


Mage Pe $(0) + 20 - (17) 

wW(B) is a relation between the maximum load per unit length on the founda- 
tion and the average deflection across the loaded width. Inasmuch as y (8) 
tends to unity as o& tends to infinity, the $" (y) integral (16) is divergent. How- 
ever, oscillating divergent infinite integrals can be considered as oscillating 
divergent series and can sometimes be summed, whereas non-oscillating 
divergent series cannot be summed. 

Add and subtract the expression / (f) in the numerator of (16): 


+ plA) i? 


Treat Integral I using Cesaro’s method as indicated by Hardy® and others: 


I 


5. “Divergent Series,” by G. H. Hardy, Oxford at the Clarendon Press, 1949, 
p. 11. 
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\ 
(14) 


20 
Given f(«) d« divergent 


Tv 
form lim (1 -~/T) d«=S 


then s= d= 


Now f(@)= cos a J and the lower limit a= 0. 
It is further stipulated that J is greater than some finite quantity which does 
not tend to zero. 


Bvaluate /(1 -~</T) cos<Jd«< by parts: 


(T sin] T Osins.0 +1 cosy T = 1 coss.0) 


(sins T-0). 


Since cos J.0=1, sin J.0= 0; and siny T and cos J T are bounded for all T. 
.. the divergent integral S= 


Note that the sumbol vy means ‘corresponds’ since the sign of equality 
implies convergence. 

Since integral I corresponds to zero, the divergent integral on the left hand 
side of (16a) is replaced by the convergent integral II. 

The same results can be shown by Abel’s method of evaluating divergent 
integrals; or, by using the exponential form indicated by Hardy.° Further- 
more, integral I in equation (16a) can also be shown to correspond to zero by 
introducing the factor e~4a@ on both sides of (16a) and allowing “a” to tend to 
zero. 

The complicated function y (g) was derived by Biot as an asymptotic 
expression for B <0.1, and tabular values were stated for B >0.1. In the 
present paper, to facilitate the required evaluation of the convergent integrals 
(14), (15), and (18) for the finite beam, w (f) is now defined as the hyperbolic 
curve 1+ 2. The value of a is determined by using the method of least 


squares,’ See Fig. 3. It can be shown that a small variation in the curve fit 
constant a has a negligible effect on the values of integrals of the types re- 
quired. A 6% variation in choosing a= 0,3 instead of 0.33, for example, re- 
sults in only a 1/2 of 1% variation in the values of integrals of the types (14), 
(15) and (18). 


6. “Divergent Series,” by G. H. Hardy, Oxford at the Clarendon Press, 1949, 
p. 12. 
7. “Graphical and Mechanical Computation,” by J. Lipka, John Wiley and 
Sons, New York, 1918, p. 126. 
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om) 
(use) 
4 
1 
foe, 


FITTING HYPERBOLIC EXPRESSION FOR TO BIOT'S TABULAR vas FOR 


-/ 

0.09 5.23 11.1 123 
Wee 142 0.08 5.68 12.5 71.0 156 5.1 
= 0.07 6.26 89.8 20h 5.7 
0.06 7.01 16.7 117 279 6.5 
0.05 8.01 20.0160 00 7.6 
0.0, 9.54 25.0 239 625 9.3 
17 

29 3h 

1.7 
Biot's/1 1.2 1.3 
Table 3 1.13 1.1 
1.04 


YA 19 asymptotic to the 
vertical axis and the 
horizontal line# (4) = 1. 


Fit the hyperbolic curve in the 
lower ranges of / to the 
logarithmic expression and the 
correspondence is good in the 
upper ranges of 4. 


1992 = 192.9 + SU97 a 
oF 0.327 say 0.33 
Method of Averages 
The? = 9 + 192.9 @ 


az 0.3 


- og 
hy 


ko 
20 
10 
— 
@M=/ 
Fig- 3 
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This simple expression for (B) permits closed form evaluations of inte- 
grals (14) and (18) when the trigonometric term vanishes. The polynomial 
expressions obtained in the denominators of the integrands in (14), (15) and 
(18), moreover, makes possible the application of the method of residues. 

The convergent integrals (14), (15) and (18) are difficult to evaluate direct- 
ly when J= (/c, }/2c is of the order of 10, due to the rapid oscillations of the 
integrands. The computational work can be considered reduced by the appli- 
cation of Euler’s transformation of series, which greatly accelerates the slow 
convergence of the required integrals. 

A more effective means of evaluation is available through the use of the 
method of residues. The values of the physical constants which are selected 
are based on representative designs of grillage beams.8,9 See Table 1 below. 


TABLE 1 
PROPERTIES OF STANDARD COLUMN GRILLAGES 


BASED ON KETCHUM'S TABLES 8 


C = 1.06 


Length I (in) (in) 


6002 2.04 7062 0.267 26 
2.68 11.6 0.232 2e2 
51-5" 508.7 2687 1368 06208 2.36 
51-10" 632-1 3.0) 14.6 0.209 2.0 
2.38 
2-0 


VI 


7'=5" 1550.3 3.56 18.7 0.191 
3010.8 4.02 22.4 0.180} 


BSRRES 


HHHHH 


= 2x 10 #/in2 for average reinf. concrete, brickwork, sand stone 


34-0" 6002 0055 
50867 0.0431) ¢ 0-490] 
5'-10" 632.1 0.0431 0.4.97 
71-5" 155043 0.0395} 
8'=11" 3010.8 0.037k} 06497 


Query 


HHHHHH 


E = 0.018 x 106 #/in2 based on sand curves 


For the standard column grillages in Ketchum’s tables the {/2c values 
appear to be fairly constant. For charts, however, the beam length is varied 
in order to illustrate the free end effects of the finite beam. As} is in- 
creased, the end effects become negligible and the beam behaves as if 
infinitely long. 1/3 


c = functional relation or fundamental length = [eca-)2)Ent 
2b = beam width 


8. “Structural Engineers’ Handbook,” by M. S. Ketchum, McGraw-Hill Book 
Company, New York, 1924, Part Il, pp. 275-283. 

9. “Progress Report of Special Committee to Codify Present Practice on 
Bearing Values of Soils for Foundations,” A.S.C.E. Papers and Discus- 
sions, Vol. 43, 1917, p. 1176. 
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Numerical values of b/c which correspond to the particular values of 
Young’s modulus E for the foundation, are introduced into the expression for 
~@). Fourth degree polynomials in o are then obtained in the denominators 
of the integrands. 

For b/c = 0.205, which corresponds to E = 2 x 106 #/in2, the oscillating 
integrals in (14), (15) and (18) are transformed into integrals containing the 
negative exponential e-JY and some additional explicit trigonometric terms. 
Numerical evaluations of damped non-oscillating integrals are much simpler 
than for oscillating integrals. 

The path chosen for contour integration is in the first quadrant along the 
real axis, quarter circle and imaginary axis. The roots of the biquadratic in 
the denominators are the singularities of the integrands. There is one in 
each quadrant and therefore only the residue of the pole in the first quadrant 


need be evaluated. 


pole at a+ bi 


x 
For the integral in (14), for example: 
2 a4 
1 
Consider 
where k,= 1.6097561, based on b/c=0.2053 ices, y(P)=1 + 
The integral along the imaginary axis is obtained by letting z = dy in (19), 


the + 4 (29) 
(y's + )2 + 


The integral along the real axis is obtained by letting z = x in (19), 


x dx + i x2 x dx (21) 
+x+ky 
The absolute value of the integral along the quarter circle is shown to tend 


to zero as R tends to infinity by letting z= R e!# in (19). Noting that 
|emto| = 1, and 


and where O< <7, 
2 


1, < forJ>0O (22) 
+R-k) 


771-10 


“4 
I; Zz 


The denominator is two degrees higher than the numerator in R and 
hence | Ig as R+oo- 


The residue=bj — p(s.) where _ g@ 
q' (zo) q\z) + 


Zo= & + bi=0-8080650 + 0.98099L5 4, 


Z=B, 


(2774) Residue (2774)(a + bi)? (a + bi) 


lta + bi)? 


Using the residue theorem, f,1@) dz = 271K, equate the real parts of I; 


and Ig to the real part of 277i ) K, where 2 K equals the sum of the residues 
at enclosed singularities. In this case, there is just one in the first quadrant. 


$l 1.6097 0.150682 sin 0480806507 ) 
+ Le 


The treatment of the integrals (15) and (18) is similar. The results are 
stated below: 


Lots -0.01810608in0.80806507) 


ely 
red 1. 1 d< 


277 98099L5S (9,2881365 cos 0.8080650F 
4022075301 sin 0.8060650 ) 


> 1. 7 ? 


For b/c = 0.044 and E = 0.018 x 10° 4 /in2 corresponding to the ‘sand’ 
foundation, the analagous identities follow: 


-R ot (27) 
+ fe + 
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(28) 


wes) = 7-5)cos d< 2772596745 (9, 3368110081 T 


@ 5 e Ty (29) 
+7 
_Jn view of the previous discussion concerning the correspondence of 


9c08 ot 7 da to zero, it is to be noted that the derivatives of (25) and (28) 
with respect to J yield (26) and (29). This serves as an additional check on 
the numerical work. 

By this procedure the oscillating integrals are summed by means of 
numerical evaluations of damped integrals and additional explicit trigonomet- 
ric terms. Thus slowly convergent oscillating integrals become rapidly con- 
vergent damped integrals. This fact arises when we consider that the larger 
the value of J the more rapidly the original functions oscillate. Rather small 
intervals in Simpson’s rule are consequently required in order to secure a 
desired accuracy. For the damped integrals, however, large values of 5 
cause the functions to decrease rapidly due to the negative exponential. 

Table 2 gives the values of the integrals in (14), (15) and (18) for those 
values of J corresponding to f/c, )/2c and 0, which are required in order to 
obtain the maximum moment by the superposition procedure. This table 
gives the results when one considers the practical designs in Ketchum,® for 
two limiting cases of elastic foundations: one of high (E = 2,000,000 # /in2) 
and one of low (E= 18,000 #/in?) Young’s modulus. 

Parametric curves of Mmax/Pc as ordinate against 9/2c as abscissa are 
shown in Fig. 5 for these two extremes of elastic foundation. The values of 
Young’s modulus are incorporated into the parameter b/c. Hence the curves 
are most general and can be applied when any of the geometric and elastic 
quantities are varied. This is shown in the section on comparison with test 
results, 

Infinitely long beams, or those sufficiently long so as to have negligible 
end effects, are represented by horizontal lines in Fig. 5. The F curves 
indicate that, other factors being equal, the moment in the beam is less for a 
foundation of large E than for one of small E. This occurs despite the fact 
that the coefficient F is somewhat greater for large than for small E. For 
4/2c = 3.4, for example, F is 0.275 as compared with 0.204. The value of c, 
however, is much larger for small E as it varies inversely with the cube root 
of E. Since this increase in c is much greater than the decrease in F, the 
maximum moment PcF is higher for the foundation of lower E. This is to be 
expected from physical considerations since the stiffer foundation should re- 
sult in smaller beam deflection and consequently less moment. 

This behavior is indicated by the graphs for all grillage beams. A practi- 
cal example of a grillage footing on sand indicates that due to the necessity 
of keeping the pressure below a certain specified value, the area of the footing 
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a 
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SAND 
Closed Forms: $(0)=0.2039169 g"(0)=0.6282867 


Eqe(1s) Eqe(15) Eq- (15) Eq. (16) 

060302 0.033136 = 002232 0.054421 0.2287 0.2075 

1.0 -0.094180 =-0.024950 =0.09876 0.029918 0.2118 

2.35 ~0.01488 0.00075516 0.02570 ~=0.00098657 0.0001526 0.2019 


5-6 9736.5 «8 0798.1 42.6 -13. 0.2039 
x10-8 x10~ x10 


ROCK 


Closed Forms: 9$(0)=0.27533 (0)=0.48237 
2635 -0.001506 0.02488 0.006709 0.2659 
0.015928 0.00038780 0.019742 0.0000927 0.2736 
34660,82 934.60) -9063.2 1016,18 -1238.566 0.2753 


x1o~ 


Note: Shear and moment are equal to zero in each case at the ends 
of the beam. 23% typical for grillage beams on rock 
in Ketchum's tables. 


Free End Effects: The maximum change in per cent of Mnygx/Pe 
Yor the finite beam as compared with the infinite beam is 
less than 


requires beam lengths such as to make {/2c > 1.19 An increase in the mo- 
ment of inertia of the beam, I, has the same effect as a decrease in the 
Young’s modulus of foundation, E, and results in an increased maximum 
moment in the beam. 

The influence of the free ends is noted for {/2c < 3.4 in the case of rock, 
and /2c < 2.35 in the case of sand foundations. The end effects are practi- 
cally the same as for an infinitely long beam. Increasing the length of the 
beam beyond /2c = 3.4 for a rock foundation does not increase the maximum 
moment, contrary to the practical design assumption that the moment varies 
directly as the beam length. This fact is indicated by the F curves and is 
often recognized intuitively in practice in designing, for example, bearing 
plates which must be extended in length for other than structural reasons. 

In this case an arbitrary equivalent length is assumed acted on by an as- 
sumed uniform bearing pressure equal to the load per unit width divided by 
this arbitrary equivalent length. 


10. “Foundations, Abutments and Footings,” by G. A. Hool and W. S. Kinne, 
McGraw-Hill Book Co., New York, 1923, pp. 245-247. 
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TABLE 2 


PARAMETRIC CURVES OF AGAINST 
BASED ON DATA IN TABLES 1 AND 2 


ROCK: 2,000,000 #/in? 
b/e = 0.205 


F=0.20) 
SAND: 18,000 #/in2, b/c = 0.0h4 


5 
Fig. 5 


Note: Point A corresponds to the designs in Ketchum's tables. 
See Table l. 


The assumption of uniform pressure gives approximately twice the maxi- 
mum moment determined from the F curves. The design of the beam is 
therefore conservative when the uniform pressure distribution is assumed. 
The maximum pressure on the foundation computed by means of the elastic 
analysis, however, is found to exceed the prescribed allowable bearing pres- 
sure considerably, as is demonstrated in the following illustrative problem. 

The elastic analysis gives considerably less flexural stress in the beam. 
Due to the higher bearing pressure under the central portion of the beam, a 


greater width of contact is desirable in order to avoid local crushing of the 
foundation. 


Maax/Pe = F 
0.50 
7 0.40 
F=0.275 
asymptotes 
0.20 
0.10 
0 
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ILLUSTRATIVE PROBLEM 


Consider a 20 I 90, 7'- 5" long, I = 1550.3 in4, b = 3.56 in., and 
S = 155.0 in3, Take Ep = 30,000 kips/in2, E= 2,000 kips/in*, Cayg = 1.06, 
) =1/4, and Ep/E = 15. Compute the fundamental length: 


es )2) [2.0602 - #) 15 2046(I/o) 
Eb 


1/3 


for all beams on rock. For the particular beam, substitute the moment of 
inertia, I, and the half flange width b. ... c = 18.6". 

Hence b/c = 3.56/18.6 = 0.191, and (/2c = 7.42 x 12/2 x 18.6 =: 2.39. Enter 
the F curve for rock in Fig. 5 with {/2c = 2.39 and obtain F = 0.265. Note: 
Most beams in Ketchum’s”® grillage tables average b/c = 0.205, and 

f/2c = 2.35. Hence there is only a slight reduction in the F value from that 
of an infinitely long beam. 


Mmax = PcF = P(18.6)(0.265) = 4.93 P. 
The units are Mma, in inch-kips, and P in kips. For a steel beam on rock 
or average concrete, Eb/E = 15 and Mma, = 0.653P z If. 
For the assumption of uniform pressure: 
P=11.1P 
4o93/11.1= 


Note: In this case, which is typical of the grillage beams in Ketchum’s 
tables. The maximum beam moment and stresses are only 44% of those 
obtained when uniform bearing pressure is assumed. 

By equating Mmax= 4.93 P, determined from the F curve, to Mmax = P/4B, 
of the spring theory, an equivalent spring constant, k, is determined. 


493 where 4 = 
A = Wieh.93=0.0507 (1/1n) 
1x30, 00021550.3x(0.0507 30,000 kips/in® 
=1230 kips/in? I= 1550.3 


Since f=7 /Bo= 62" is less than the beam length, the tendency towards uplift 
is assumed prevented by the foundation sustaining tension. Using the spring 


3 -De- 
theory, shear V= = *cosA,x, 
P_av _P Ax (sin +cos Ax)= ky. 


The function = Rx (sin bx + cos fox ) has been tabulated by Timoshenko! ! 

and Hetenyi.!2 

11, “Strength of Materials,” by S. Timoshenko, D. Van Nostrand Company, 
New York, 1941, Part Il, p. 5. 

12. “Beams on Elastic Foundation,” by M. Hetenyi, U. of Michigan Press, 
Ann Arbor, Michigan, 1946, p. 217. 
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Now, under the load, x = 0 and p= pmax = P 8, /2 = P(0.0507)/2 
=0.0254 P kips per inch of beam length. 
”. max unit pressure = Pmax/2b= 0.0254P/7.13 = 0.00356P 
(kips/in?) 
For a maximum allowable pressure of 0.500 kips/in2, the maximum load 
P= 140 kips. 
For a maximum allowable bending stress of 18 kips/in?, M=fS or 


4.93P = 18 x 155.0, P= 566 kips 
Using this elastic analysis, the allowable bearing pressure governs the 
maximum load to be imposed on the beam. Under the uniform pressure 


assumption 


Pallow.brg.= 0.500x7 «13x(7 o42x12)= 317 kips 


Pallow.bend.= x8 — 251 kips 


Here the maximum bending stress in the beam governs the maximum load 


to be imposed. 
Now at Ax =1, however, the maximum pressure drops to 50.8% of its 
value under the load.!1 x= 1/0.0507 = 19.7". Lower values occur beyond a 


distance equal to about one fourth the beam length from the load. 

In this example it is shown that if the beam is selected based on a specified 
permissible uniform bearing pressure, the flexural stress computed by the 
elastic analysis is low (about 1/2), and the actual maximum bearing pressure 
is higher (about twice) than the original permissible bearing pressure. The 
deflection under the load using the spring theory is: 


0.0507P 
P B,/2k 2x 1230" 0.0000206P inches, where P is in kips. 


(P/2)(f/2)8 


Under the uniform bearing assumption, max. defl. = BEI 


7.423 x 1728P 
max, defl, = @ x 16 x 30,000 x 1550.3 = 0.000119P inches. 


Here less deflection (about 1/6) is indicated when the spring theory is used. 
The particular beam selected is adequate to sustain a concentrated load of 
140 kips when the elastic analysis is used. The assumption of uniform pres- 
sure indicates a permissible load of 250 kips. The results show that from 
purely elastic considerations the usual treatment of the grillage beam, the 
primary function of which is to reduce the maximum pressure on the founda- 
tion, is considerably in error. 


Comparison with Test Results 


It is interesting to compare the maximum stress determined by means of 
the curves of Fig. 5 with some test results stated by de Beer.13,14 The 


13, “Calcul de Poutres reposent sur le sol,” by E. E. de Beer, Ann. Trav. 
publics Belg., Vol. 49, 1948, pp. 393-413, 525-552, 653-691, 721-749, 
especially p. 740 and p. 746. 

14. “Tests for the Determination of the Distribution of Soil Reactions Under- 
neath Beams Resting on Soil,” by E. de Beer, Proc. Sec. Intern. Conf. 

Soil Mechanics and Found, Eng., Vol. 2, 1948, pp. 142-148. 
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Boussinesq solution was employed by de Beer to compute the stresses in the 
foundation for the case of an ideal soil treated as homogeneous and isotropic. 
De Beer used a method of successive approximations to secure a nearly coin- 
cident beam and soil displacement pattern. Since de Beer’s results are a 
function of the parameter width/length of the loaded area of the foundation, 
the results obtained are not general but must be repeated for each ratio. 

The parametric curves of the present paper were originally computed to 
encompass the grillage beam designs shown by Ketchum.® The choices of 
the parameter b/c were originally for the cases of steel beams on rock and 
for steel beams on ideal sand. The parametric form of these curves, how- 
ever, allows ready substitution of de Beer’s numerical values for the physi- 
cal quantities concerned. The maximum moment can be quickly estimated. 

Use de Beer’s values for a stiff T beam: Ixx = 949 cm4; Section 
Modulus = 353 cm3; Ey = 2 x 108 kg/cm?2; /=1 meter; width = 0.25 m.; and 
Esand = 544 kg/cm? (experimental average). Enter the parametric curves of 
the present paper, Fig. 5, with: 


e= [on & (2 3°) | V3 


c= 60.8 em; 100 0.825 
2260.8 


It is seen by linear extrapolation of the rock curve that F = Mmax/Pc = 0.25 

and Mmax = PcF = 60.8 x 0.25P = 15.2P cm-kg. For P = 2500 kg, which is an 

average load value in de Beer’s tables, 0=15.2 x 2500/353 = 108 kg/cm2. 
stated by de Beer 


Stiff Tee Unif. Pressure This ‘K Method 2nd Degree Strain 


Beam Assumption Paper Parebola Gages 
stress 88 5 108 101 112.2 
kg/cm 


In the case of the flexible T beam considered by de Beer the moment of 
inertia is reduced to 54.7 cm*, section modulus = 48.7 cm® and the other 
quantities remain the same. The results for P = 2500 kg are as follows: 


stated by de Beer 


Flexible Unif.Pressure This ‘K Method 2nd Degree Strain 
Tee Bean Assumtion Paper Parabola Gages 


CONC LUSIONS 


The parametric curves of Mmax/Pc against {/2c shown are for two values 
of b/c which correspond to the numerical values of Erock = 2,000,000 #/in* 
and Egand = 18,000 #/in“. From these curves determinations are made of 
the maximum beam moments which do not contain the spring constant. For 
any significant change in one or more of the elastic or geometric quantities 
embodied in the ratios M/Pc, §/2c and b/c, the parametric curves can be 
used, It is noted that for the entire range of grillage beam designs shown in 
Ketchum’s tables ,® the value of b/c is substantially constant for any 
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particular value of E. Using E = 2,000,000 #/in2, for example, b/c = 0.205. 

A non-experimental determination of k is readily obtained by equating the 
maximum moment in the spring foundation. Using this value of k, basic 
problems, such as grillage beams, can be treated and their solution compared 
with results obtained by other means. Instead of using expressions for de- 
flection, slope, moment and shear containing an experimental constant k, the 
formulas of this paper employ the elastic constants E and ) which are experi- 
mentally determined. The known values of E and ) characterize only the 
foundation material and are independent of the structure to be imposed. 

By use of the expression for the moment in the infinite beam under a 
single concentrated load the maximum moment for a beam of arbitrary length 
subjected to any loading can be developed by superposition. 

The computations of this paper indicate that the greater pressures present 
in the vicinity of the concentration result in considerably less bending stress 
in the beam. The larger maximum foundation pressure is not substantially 
reduced by the prevalent practice of increasing the beam length. 


APPENDIX. NOTATION 

a (dimensionless) = Ac 

b (in.)= 1/2 flange width of beam 

b/c (dimensionless) = parameter of problem: ratio of 1/2 flange width of 
beam, 2b, to afundamental length, c. This parameter is substan- 
tially constant for typical grillage tables for each type of foundation 
(i.e., for any particular E). 

B (dimensionless) = bo/c parameter. 


Bo (1/in.) k/4E I reciprocal of the fundamental length arising in the 
spring theory. Sometimes termed ‘characteristic’ or ‘damping 


factor.’ 
c (in.) = fundamental length (1 -9 2) 


C(dimensionless) = coefficient representing effect of changing distribution of 
loading across width 2b. C =1.13 for uniform deflection across 
width 2b. C =1.06 used as average in the computations. C =1 for 
uniform pressure across width 2b. 

Ep (#/in2)= Young’s modulus for the beam. 

E (#/in2) = Young’s modulus for the foundation. 

F = Mmax/Pc = dimensionless quantity by which Pc is multiplied to obtain, 
Mmax, the moment under the concentration P. F is determined 
after entering the curves in Fig. 5 with values of §/2c and b/c. 

f (#/in2) = allowable stress. 

f= x/e = parameter 

= moment of inertia of the beam. 

k (#/in. length/ in. deflection) = spring modulus. 

ky @#/in3: 4 /in. length/in. deflection/unit width of beam) = spring modulus of 
the foundation based on unit width of beam: sometimes termed 
‘coefficient of subgrade reaction.’ 

§ (in.) = length of finite beam. 

Lmax (#/in. of beam length) = 2q,> = maximum load on the foundation surface. 

Lavg. (#/in. of beam length) = average load on the foundation surface. 

Q/c, §/2c = parameter of problem. 
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A (1/in.) = reciprocal of fundamental wave length of loading in beam length 
direction. 

M (in-¥ ) = moment in beam due to load P. 

M,(in-#) = unknown annulling moments applied at those points of the infinite 
beam, corresponding to the end points of the finite beam, so as 
to give the free end conditions. 

) (dimensionless) = Poisson’s ratio of the foundation. This is defined as the 
ratio of lateral contraction to longitudinal extension for a 
particular material. Assumed constant = 1/4 in computations. 

P (#) = concentrated load at center of beam. 

p (#/in. of beam length) = reaction of beam on foundation. 

p(x) (#/in. of beam length) = load distribution along the beam. 

Qo (#) = unknown annulling shear forces applied at those points of the infinite 
beam, corresponding to the end points of the finite beam, so 
as to give the free end conditions. 

do (4/in2) = prescribed constant pressure distribution across the flange width. 

S (in3) = section modulus of beam. 

V (#) = dM/dx = shear. 

Wo(y) (in.) = maximum deflection of foundation surface along a section across 
the flange width of beam. 

Wavg. (in.) = average deflection of the foundation surface across the flange. 


ad = dimensionless quantity by which Pc 
(f= is multiplied to obtain the moment at 
any point in an infinite beam due to a 

concentrated load, P, at the origin. 


dimensionless quantity by which P is 
At multiplied to obtain the shear at any 


ax point in an infinite beam due to a con- 
centrated load, P, at the origin; 
dimensionless quantity by which the 
moment M, is multiplied to obtain the 
moment at any point in an infinite 
beam due to Mo at the origin. 


2 
2") = dimensionless quantity by which M,/c 
is multiplied to obtain the shear at 
any point in an infinite beam due to 
Mj at the origin. 


x (in.) = distance of any point along the beam from the origin. 

y (in.) = vertical deflection of any point along the beam. 

y (8) (dimensionless) = a hyperbolic function obtained by fitting by means of 
least squares to Biot’s asymptotic expression and tabular 


values. It relates Lmax to Wavg 


| 
| 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are presented below. Technical-division 
sponsorship is indicated by an abbreviation at the end of each Paper Number, the symbols 
referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics 
(EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), Sanitary 
Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping 
(SU), and Waterways (WW) divisions. For titles and order coupons, refer to the appropriate 
issue of “Civil Engineering” or write for a cumulative price list. 


VOLUME 80 (1954) 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 
475(SM), 476(SM), 477(SM), 478(SM)°, 479(HY)°, 480(ST)°, 481(SA)°, 482(HY), 483(HY). 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)°, 488(ST)°, 489(HY), 490(HY), 491(HY)°, 
492(SA), 493(SA), 494(SA), 495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 500(HW), 501(HW)°, 
503(WW), 504(Ww)°, 505(CO), 506(CO)*, 507(CP), 508(CP), 509(CP), 510(CP), 

CP). 


OCTOBER: 512(SM), 513(SM), 514(SM), 515(SM), 516(SM), 517(PO), 518(SM)°, 519(IR), 520(R), 
521(IR), 522(IR)°, 523(AT)°, 524(SU), 525(SU)°, 526(EM), 527(EM), 528(EM), 529(EM), 
530(EM)°, 531(EM), 532(EM)~, 533(PO). 


NOVEMBER: 534(HY), 535(HY), 536(HY), 537(HY), 538(HY), 539(ST), 540(ST), 541(ST), 542(ST), 
543(ST), 544(ST), 545(SA), 546(SA), 547(SA), 548(SM), 549(SM), 550(SM), 551(SM), 552(SA), 
553(SM)°, 554(SA), 555(SA), 556(SA), 557(SA). 


DECEMBER: 558(ST), 559(ST), 560(ST), 561(ST), 562(ST), 563(ST)°, 564(HY), 565(HY), 566(HY), 
567(HY), 568(HY)°, 569(SM), 570(SM), 571(SM), 572(SM)°, 573(SM)°, 574(SU), 575(SU), 576(SU), 
577(SU), 578(HY), 579(ST), 580(SU), 581(SU), 582(Index). 


VOLUME 81 (1955) 


JANUARY: 583(ST), 584(ST), 585(ST), 586(ST), 587(ST), 588(ST), 589(ST)°, 590(SA), 591(SA), 
592(SA), 593(SA), 594(SA), 595(SA)°, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(WW), 
616(WW), 617(IR), 618(IR), 619(1R), 620(1R), 621(IR)©, 622(IR), 623(1R), 624(HY)®, 625(HY), 
626(HY), 627(HY), 628(HY), 629(HY), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(PO), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)°, 642(SA), 
643(SA), 644(SA), 645(SA), 646(SA), 647(SA)°, 648(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)©, 655(SA), 656(SM)°, 657(SM)°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)°, 662(ST), 663(ST), 664(ST)°, 665(HY)°, 666(HY), 667(HY), 
ceaty), 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 


MAY: 679(ST), 680(ST), 681(ST), 682(ST)°, 683(ST), 684(ST), 685(SA), 686(SA), 687(SA), 688(SA), 
689(SA)°, 690(EM), 691(EM), 692(EM), 693(EM), 694(EM), 695(EM), 696(PO), 697(PO), 698(SA), 
699(PO)©, 700(PO), 791(ST)°. 


JUNE: 702(HW), 703(HW), 704(HW)°, 705(IR), 706(IR), 707(IR), 708(IR), 709(HY)°, 710(CP), 
T11(CP), 712(CP), 713(CP)©, 714(HY), 715(HY), 716(HY), 717(HY), 718(SM)°, 719(HY)¢, 
720(AT), 721(AT), 722(SU), 723(WW), 724(WW), 725(WW), 726(WW)°, 727(WW), 728(IR), 
729(IR), 730(SU)°, 731(SU). 


JULY: 732(ST), 733(ST), 734(ST), 735(ST), 736(ST), 737(PO), 738(PO), 739(PO), 740(PO), 
741(PO), 742(PO), 743(HY), 744(HY), 745(HY), 746(HY), 747(HY), 748(HY)°, 749(SA), 750(SA), 
751(SA), 752(SA)°, 753(SM), 754(SM), 755(SM), 756(SM), 757(SM), 758(CO)©, 759(SM)°, 
760(ww)”. 


AUGUST: 761(BD), 762(ST), 763(ST), 764(ST), 765(ST)©, 766(CP), 767(CP), 768(CP), 769(CP), 
770(CP), 771(EM), 772(EM), 773(SA), 774(EM), 775(EM), 776(EM)©, 777(AT), 778(AT), 
779(SA), 780(SA), 781(SA), 782(SA)°, 783(HW), 784(HW), 785(CP), 786(ST). 


c. Discussion of several papers, grouped by Divisions. 
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